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ABSTRACT 

We use the results of the SCUBA Local Universe Galaxy Survey, a submillimetre survey 
of galaxies in the nearby Universe, to investigate the relationship between the far-infrared- 
submillimetre and radio emission of galaxies at both low and high redshift. At low redshift 
we show that the correlation between radio and far-infrared emission is much stronger than 
the correlation between radio and submillimetre emission, which is evidence that massive 
stars are the source of both the far-infrared and radio emission. At high redshift we show that 
the submillimetre sources detected by SCUBA are brighter sources of radio emission than are 
predicted from the properties of galaxies in the local Universe. We discuss possible reasons 
for the cosmic evolution of the relationship between radio and far-infrared emission. 



Key words: 

galaxies. 



surveys - galaxies: ISM - galaxies: distances and redshifts - radio continuum: 



1 INTRODUCTION 

One of the most intriguing relationships in astronomy is that be- 
tween the non-thermal radio and far-infrared (FIR) emission from 
galaxies. It is one of the strongest correlations in astronomy, with 
a tight correlation between FIR and radio luminosity over five 
decades of luminosity (e.g. Price & Duric 1992; Yun, Reddy & 
Condon 2001). The relationship holds for a remarkably wide va- 
riety of galaxy types, and is observed not only globally but also 
on ~kiloparsec scales within the disks of individual galaxies (e.g. 
Boulanger & Perault 1988; Bicay & Helou 1990; Fitt et al. 1992; 
Murphy et al. 2006). Some authors also find evidence of the rela- 
tionship at high redshifts (e.g. Garrett 2002; Appleton et al. 2004; 
Kovacs et al. 2006). Yet, although the relationship has been known 
for over two decades (e.g. Condon 1992 and references therein), 
the cause of the relationship is still unclear. 

The standard explanation of the relationship is that both the 
FIR and radio emission are caused by high-mass stars, with the 
stars both heating the dust, which then emits FIR emission, and 
producing the relativistic electrons responsible for synchrotron ra- 
diation when at the end of their lives they explode as supernovae 
(see e.g. Condon 1992 and references therein). Alternative explana- 
tions have been proposed. However, most workers in the field have 
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assumed that the standard explanation (i.e. ongoing star-formation) 
is the correct one. 

Nevertheless, despite its widespread acceptance, there is sur- 
prisingly little direct evidence for the standard explanation and it 
has one major problem. The problem is that the slope of the FIR- 
radio relation is roughly unity (Helou, Soifer & Rowan-Robinson 
1985; Devereux & Bales 1989) whereas the synchrotron emissivity 
does not just depend on the number-density of relativistic electrons 
but also on the strength of the magnetic field. Although there are 
only a few measurements from Zeeman splitting of the magnetic 
fields in galaxies (Thompson et al. 2006 and references therein), 
it is possible to estimate the magnetic field in a galaxy using the 
'minimum energy' argument (Longair 1983), and estimates of the 
magnetic field made in this way can vary by as much as a factor of 
100 from normal galaxies like our own to vigorously star-forming 
galaxies such as Arp 220 (Thompson et al. 2006). The radio spec- 
trum of a galaxy roughly follows a power-law (5* oc and the 
synchrotron emissivity is proportional to B^""''". The radio emis- 
sivity in Arp 220 will therefore be much higher than in our own 
galaxy, not only because of the much greater production rate of 
relativistic electrons per cubic kiloparsec, but also because of the 
stronger magnetic field. However, it is still possible to explain the 
unity slope of the FIR-radio relation if the relativistic electrons lose 
all their energy within their galaxy, and thus an electron in Arp 220 
loses its energy at a faster rate but ultimately only loses the same 
amount of energy as in a galaxy like our own. 

This scenario is described by the 'calorimeter' model (Volk 
1989; Lisenfeld, Volk & Xu 1996), in which the dust-heating stel- 
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lar UV photons and the relativistic electrons lose most of their en- 
ergy within the galaxy. The calorimeter model will only give a unity 
slope for the FIR-radio relation if the ratio of inverse Compton to 
synchrotron losses is roughly the same in all galaxies, which leads 
to the prediction that the energy density of the interstellar radia- 
tion field and that of the magnetic field have a constant ratio for 
galaxies. Lisenfeld et al. (1996) used observations of 114 'normal' 
late-type galaxies to estimate the mean energy densities of the inter- 
stellar radiation field and, using the 'minimum energy' argument, 
the magnetic field. They showed a linear correlation between the 
two, suggesting turbulent dynamo activity in galaxies as a possible 
explanation. Calorimeter-type models have also been presented by 
several other authors (e.g. Pohl 1994). 

A more fundamental problem with the calorimeter model is 
that the radio spectra of star-forming galaxies do not exhibit the 
steepening that one would expect if the electrons lose all their en- 
ergy within the galaxy in which they were formed. A recent inter- 
esting suggestion has been that the magnetic fields in galaxies are 
actually much higher than the estimates from the minimum energy 
argument, which would definitely lead to an electron losing all its 
energy before it diffuses out of its galaxy (Thompson et al. 2006). 
Thompson et al. explain the lack of clear curvature in the radio 
spectra as the result of ionization losses, which flatten the spectrum 
at low frequencies. 

An alternative idea is that the FIR-radio relation holds for the 
'optically thin' scenario in which the majority of the stellar UV 
photons and relativistic electrons do escape from the galaxy disk 
(e.g. Helou & Bicay 1993). The FIR-radio correlation is then at- 
tributed to a constant ratio between the production rates of dust- 
heating UV photons and synchrotron-emitting relativistic electrons. 

Niklas & Beck (1997) argue that observations do not support 
either the calorimeter or 'optically-thin' model and suggest an al- 
ternative model in which the FIR-radio relation arises from a com- 
bination of correlations; that between the volume density of cool 
(neutral) gas and the star formation rate and magnetic field strength, 
and that between the star formation rate and the FIR luminosity, as- 
suming equipartion between the cosmic ray energy density and the 
magnetic field strength. In this scenario the cool neutral gas clouds 
play a central role. Such a model, unlike either the calorimeter (op- 
tically thick) or optically thin models, can explain the fact that the 
FIR-radio relation is not only observed globally but is also ob- 
served on kiloparsec scales within the disks of individual galaxies. 

Alternative explanations have also been proposed by a num- 
ber of other authors. Bettens et al. (1993) proposed that the cosmic 
ray flux, which as in the 'standard explanation' produces the non- 
thermal radio emission, drives the evolution of molecular clouds; a 
higher flux of cosmic rays produces a higher ionization fraction in 
the cloud, which delays the cloud's collapse, and ultimately leads 
to a larger fraction of OB stars and more FIR emission from the sur- 
rounding dust. However, other authors (e.g. Niklas & Beck 1997, 
Murgia et al. 2005) have pointed out problems with this model; for 
example it fails to take into account the effect of the the interstellar 
magnetic field strength on the synchrotron luminosity. 

Murgia et al. (2005) proposed a model in which the relation 
between CO, radio continuum and FIR emission results from hy- 
drostatic pressure in galaxies, which can explain the correlation on 
both global and local scales. In this model the relationship between 
the CO and radio is explained by hydrostatic pressure acting as a 
regulating mechanism between the molecular gas fraction and the 
magnetic field strength, while the link to the FIR emission comes 
from models by Dopita et al. (2004) which connect FIR emission 



with interstellar pressure. This model thus avoids an explicit depen- 
dence on star formation. 

Most of the various models proposed to explain the FIR-radio 
relationship require a link between the gas density and the magnetic 
field strength. One suggestion as to possible physical processes be- 
hind such a model is magnetohydrodynamic turbulence (Groves et 
al. 2003). 

Whatever its cause, the strong correlation between FIR and 
radio emission in nearby galaxies provided submillimetre (submm) 
astronomers with a valuable tool for estimating the redshifts of the 
galaxies discovered in the first deep submm surveys (Small, Ivison 
& Blain 1997; Hughes et al. 1998; Eales et al. 1999). These galax- 
ies are very faint at optical wavelengths, which has made it very 
difficult to measure redshifts using conventional spectroscopy, and 
until recently very few spectroscopic redshifts existed (Chapman et 
al. 2005). Carilli & Yun (1999) first suggested that an alternative to 
optical spectroscopy would be to assume that the FIR-radio rela- 
tion is the same at all redshifts, and then estimate redshifts from the 
ratio of submm to radio flux densities. 

In this paper we explore the FIR-radio relationship at both 
high and low redshifts, using the results of the recently-completed 
SCUBA Local Universe Galaxy Survey (SLUGS). This is a submm 
survey of ~ 200 galaxies in the local Universe, half drawn from a 
sample selected in the FIR (Dunne et al. 2000b; herein the IRS sam- 
ple) and half drawn from an optically-selected (herein OS) sample 
(Vlahakis, Dunne & Eales 2005). Our submm measurements allow 
us to test a basic prediction of the standard theory. Because of the 
greater intensity of the interstellar radiation field in regions con- 
taining large numbers of OB stars, the dust in these regions will 
be hotter than the dust in the general interstellar medium. There- 
fore, the standard hypothesis makes the prediction that the cor- 
relation between FIR and radio emission will be tighter than the 
correlation between submm and radio emission. Previous authors 
have found evidence both in favour of and against this prediction 
by decomposing the FIR emission into warm and cold components 
(e.g. Hoemes, Berkhuijsen and Xu 1998; Pierini et al. 2003), but 
this previous work was based on FIR measurements alone rather 
than measurements at submm wavelengths. In Section [2] we use 
our SLUGS data to test the prediction. In Section [3] we then con- 
sider the FIR-radio relationship at high redshift. Because SLUGS 
is the only large sample of galaxies for which there are measured 
spectral energy distributions extending through the FIR and submm 
wavebands, we can predict the ratios of radio to submm emission 
one would see at high redshift ;/ there is no cosmic evolution in 
the FIR-radio relationship. Dunne, Clements & Eales (2000) per- 
formed a similar analysis using the IRS data alone. Using the spec- 
troscopic redshifts from Chapman et al. (2005), we show that the 
sources detected in deep SCUBA surveys are generally brighter ra- 
dio sources than one would predict from the FIR-radio relationship 
seen in the local Universe. 



2 THE FIR-RADIO RELATIONSHIP IN THE 
LOW-REDSHIFT UNIVERSE 

The radio parameters for the OS sample are given in Table [T] Ra- 
dio (1.4 GHz) fluxes were taken from the literature, most from the 
NRAO/VLA Sky Survey (NVSS; Condon et al. 1998) or otherwise 
from the VLA FIRST (Faint Images of the Radio Sky at Twenty 
Centimeters) survey (Becker, White & Helfand 1995). For objects 
without literature radio fluxes the NVSS radio images were used 
to either measure the flux of a faint source or, where no source was 
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Table 1. Literature 1.4 GHz Radio fluxes. Table 1 - continued 



(1) 


(2) 


(3) 


(4) 


(1) (2) (3) 


(4) 


Name 


51.4 

(mJy) 


log L\A 
(WHz-lsr-i) 




Name 5i.4 logLi.4 
(mJy) (WHz-isr-i) 


„850 
"1.4 


UGC 148 


17.8 


20.71 


0.21 


UGC 10486 <1.5 < 19.96 




NGC 99 


8.2 


20.58 


0.37 


NGC 6196 1.0 20.17 


<0.57 


PGC 3563 


4.1 


20.31 


0.34 


UGC 10500 <1.4 < 19.79 




NGC 786 


6.1 


20.31 


0.43 


IC5090 32.6 21.67 


0.23 


NGC 803 


<1.5 


< 19.02 


>0.75 


IC 1368 23.7 20.80 


0.12 


UGC 5129 


1.8 


19.68 


<0.54 


NGC 7047 3.9 20.31 


0.48 


NGC 2954 


<0.84 


<19.30 




NGC 7081 18.2 20.50 


0.16 


UGC 5342 


3.8 


20.11 


0.39 


NGC 7280 <1.2 <18.82 




PGC 29536 


1.0 


20.15 


<0.68 


NGC 7442 7.8 20.83 


0.32 


NGC 3209 


1.11 


19.84 


<0.54 


NGC 7448 81.5 20.80 


0.16 


NGC 3270 


6.2 


20.60 


0.41 


NGC 7461 <1.4 <19.62 




NGC 3323 


13.4 


20.76 


0.30 


NGC 7463 16.3 20.16 




NGC 3689 


27.1 


20.52 


0.24 


111 ZW 093 5.3 21.30 


<0.29 


UGC 6496 


1.5 


19.98 


0.47 


111 ZW 095 1.5 20.14 


<0.46 


PGC 35952 


3.8 


19.99 


0.47 


UGC 12519 4.6 20.16 


0.51 


NGC 3799/3800 


49.3 


20.94 


0.18 


NGC 7653 92.2 21.43 


0.04 


NGC 3812 


1.96 


19.62 


<0.54 


NGC 7691 4.5 20.08 


<0.31 


NGC 3815 


3.2 


19.85 


0.47 


NGC 7711 <1.2 <19.51 




NGC 3920 


7.6 


20.21 


0.27 


NGC 7722 4.7 20.09 


0.47 


NGC 3987 


56.3 


21.27 


0.22 






NGC 3997 


5.2 


20.28 


<0.27 






NGC 4005 


2.4 


19.89 


<0.33 


(1) Galaxy name; (2) 1.4 GHz radio flux taken from the lit- 


NGC 4015 


<1.5 


< 19.66 




erature; NVSS (most) or FIRST; (3) 1.4 GHz luminosity; (4) 


UGC 7115 


87.3 


21.82 


—0.10 


radio-submm spectral index (a^^^), calculated as described in 


UuC 7157 


<1.4 


< 19.92 




Section|3] 




iL- /y/ 


i i.j 


1 Q 01 


U.30 








1 4 


1 Q no 








NGC 4712 


<1.2 


< 19.57 


>0.81 


seen, to measure an upper limit. The 1 .4 GHz luminosities were cal- 


PGC 47122 


<1.8 


<20.17 




culated using the standard equation (e.g. Condon et al. 1990). Apart 


MRK 1365 


23.0 


21.06 


0.06 


from a few measured radio spectral indices we usually assumed a 


UGC 8872 


<1.1 


< 19.74 




value of 0.7 (Condon 1992). For the IRS sample, radio parame- 


UGC 8883 


3.0 


20.18 


<0.47 


ters are those used in Dunne et al. (2000a), which were taken from 


UGC 8902 


12.4 


21.08 


0.31 


Condon et al. (1990); where there were flux measurements in the 


IC 979 


<1.5 


<20.17 


>0.66 


literature at more than one wavelength spectral indices were calcu- 


UGC 9110 
NGC 5522 


7.3 
17.6 


20.41 
20.78 

ZU. /D 


<0.34 
0.26 

n on 
u.zu 


lated, and otherwise a value of 0.7 was assumed. 


The relationship between 1.4 GHz luminosity and 


NGC 5980 


25 3 


20 84 


0.42 


60 (xm luminosity is shown in Figure [T] The filled points in- 


IC 1 174 


<1.5 


< 19.73 


>0.51 


dicate OS objects detected at 850 \x.vs\ while for comparison the 


UGC 10200 


11.2 


19.85 


<0.11 


open points represent the OS non-detections at 850 i^m (see Vla- 


UGC 10205 


2.9 


20.31 


0.55 


hakis et al. 2005). The IRS galaxies are plotted as cross symbols. 


NGC 6090 


48.0 


21.79 


0.12 


In the following analysis only OS galaxies with detections at the 


NGC 6103 


5.9 


20.94 


0.40 


relevant wavelengths and not upper limits 


are considered. For 


NGC 6104 


6.4 


20.87 


<0.30 


the OS sample there is a very tight correlation between 1.4 GHz 


IC 1211 


<1.3 


<19.83 


>0.56 


luminosity and 60 nm luminosity (correlation coefficient rs= 0.97, 


UGC 10325 


1 1.2 


20.77 


0.24 


probability<0.001), with a less tight correlation for the IRS 


NGC 6127 

OiZO 


<1.2 

33. / 


< 19.66 

Zl.D/ 


>0.78 
n 1 T 


sample. Figure [T] also shows the linear least- 


squares fit to the OS 


sample alone (solid line), where the slope of the log-log plot is 


NGC 6131 


6.5 


20.44 


0.39 


1.04±0.04. We note that the slope for the OS sample, both for 


NGC 6137 


438.7 


22.80 


-0.50 


Figure [T] and Figure |2l was calculated usinj 


; only those galaxies 


NGC 6146 


162.7 


22.32 


-0.32 


with detections and not upper limits, and thus this may affect the 


NGC 6154 


1.2 


19.85 


<0.64 


exact value of the slope. Apart from this, any apparent discrepancy 


NGC 6155 


13.9 


20.12 


0.08 


between our value of the slope and that found by previous authors 


UGC 10407 


18.4 


21.33 


0.06 


(e.g. Cox et al. 1988 and Devereux & Eales 1989) could also 


NGC 6166 


3713.0 


23.70 


-0.72 


be due to the fact that a) we plot 60 \x.m luminosity whereas the 


NGC 6173 


6.9 


20.94 


<0.23 


previous authors used the FIR luminosity, and b) our SLUGS 


NGC 6189 
NGC 6190 


7.9 
4.8 


20.61 
19.94 


0.40 
0.55 


sample includes higher 60 \ssa. and FIR luminosities than these 


NGC 6185 


62.1 


22.04 


< -0.13 


previous studies. 
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Figure 1.1.4 GHz luminosity versus 60 [am luminosity for the OS and IRS 
samples (circles and crosses respectively; for the OS galaxies filled cir- 
cles indicate those detected at 850 ytm while open circles indicate the non- 
detections at 850 ytm (see Vlahakis et al. (2005)); for clarity, OS galaxies 
with both 60 \xm and radio upper limits are not shown). The solid line shows 
a linear least-squares fit to the OS galaxies only, where the slope of the 
log-log plot is 1.04±0.04. The 7 OS outliers which have 'excess' 1.4 GHz 
emission relative to this correlation are all early-type galaxies. 



The detected OS galaxies appear to follow an even tighter 
60 fi.m~-radio relation than the IRS galaxies. The reason for this 
is not clear, but it could be related to the fact that the luminous 
IRAS galaxies tend to be mergers or more active systems. Another 
possible explanation could be that the several objects in the IRS 
sample which appear to be under-luminous in the radio (including 
NGC4418; see discussion by Roussel et al. 2003) are undergoing 
a very recent burst of intense star-formation and so there has not 
been time for the number of relativistic electrons to reach the equi- 
librium value. We note, however, that the apparent difference in 
scatter may also be explained by the large number of upper limits 
for the OS galaxies. 

There are 7 extreme outliers (all OS galaxies) which have 
excess 1 .4 GHz emission relative to this correlation, all of which 
are early-types (the furthest outlier is NGC6166). The 'excess' 
1.4 GHz emission most likely indicates the presence of an AGN 
in these galaxies. 

The relationship between 1.4 GHz luminosity and 
850 \im luminosity is shown in Figure |2] For the OS sample 
there is a correlation (rs=0.73, probability<0.001), but there 
is clearly a much tighter correlation between radio and 60- vim 
luminosity than between radio and 850- [xm luminosity (the rms 
scatter is 0.13 for 60 |j.m-radio and 0.40 for 850 (xm-radio). This 
is exactly the behaviour we would expect if the standard model 
is correct (see Section [TJ. We would expect to see a stronger 
correlation between the 60-iim and radio luminosity, because 
the 60-iim emission is dominated by emission from hot dust in 
star-formation regions mainly heated by OB stars, whereas the 
850-ii.m emission traces colder dust, which is heated by the general 
interstellar radiation field, which must include a component from 
older stellar populations. Vlahakis et al. (2005) showed that the 
OS galaxies contain large proportions of cold dust, with up to 
three orders of magnitude more mass of cold dust than warm dust. 
The fact that a much larger scatter at 850 \im is precisely what 
we see thus provides evidence that the standard explanation for 
the FIR~radio luminosity correlation is correct and that it is the 




10^° 10" 10'^ 10" 

L,5„ {W Hz-'sr-') 

Figure 2. 1 .4 GHz luminosity versus 850 |am luminosity for the OS and 
IRS samples (filled circles and crosses respectively; for clarity, OS galaxies 
with both submm and radio upper limits are not shown). The solid line 
shows a linear fit to the OS data alone, where the slope of the log-log plot 
is 1.08±0.16. There is a correlation between radio and 850 |am luminosity 
for the OS sample (rs =0.71, probability<0.001) but there is much more 
scatter than for the relation with 60 \xm luminosity. 

population of young (OB) stars which both heats the dust and 
provides the source of radio emission. 

Comparing the OS galaxies to the IRS galaxies there are some 
interesting points to note from Figure (2] i) Both the OS and IRS 
samples follow a similar slope (the line in Figure [2] shows the lin- 
ear fit for the OS sample alone, where the slope of the log-log plot 
is 1.08±0.16), but while for the OS sample the scatter is much 
greater than seen for the 60ii-m-radio relation (with implications 
which are discussed above), for the IRS sample there is a similar 
amount of scatter for both 60-ii.m and 850- nm luminosity with the 
radio (~0.2). ii) There is a much larger scatter for the OS sam- 
ple than the IRS sample (0.22 for the IRS sample and 0.40 for the 
OS sample, excluding those galaxies with 'excess' radio emission 
which are likely AGN-related). We can explain both these facts if 
for the IRS galaxies warm dust dominates the emission at both 60- 
and 850-nm (see Vlahakis et al. 2005). 



3 THE FIR-RADIO RELATIONSHIP IN THE 
HIGH-REDSHIFT UNIVERSE 

Carilli & Yun (1999) proposed a method of using the redshift- 
sensitive nature of the radio-submm flux density ratio as a redshift 
estimator. The flux density ratio S'85o/<S'i.4 is related to the radio- 
submm spectral index between 1.4 GHz and 850 yuca {oi^2) by the 
equation 

a?^° = 0.42 X log (^l^i) (1) 

where Sgso and 51.4 are the flux densities at 850 \im and 1.4 GHz 
respectively (as defined in Carilli & Yun 1999). Carilli & Yun 
(2000) predicted the evolution of ci^2 with redshift for 17 low- 
redshift star-forming galaxies by making polynomial fits to the ob- 
served FIR-radio SEDs of these galaxies. Likewise, Dunne et al. 
(2000a) predicted the change of with redshift for the 104 
galaxies in the IRS SLUGS sample. This method is based on the 
assumption that the FIR-radio luminosity relation is the same at 
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low and high redshifts, i.e. that more distant galaxies have simi- 
lar properties (e.g. magnetic field strength or dust grain properties) 
to those observed in the local Universe. As discussed in Carilli & 
Yun (2000), other factors which could potentially affect the relia- 
bility of this method include the spectral index assumed for the ra- 
dio synchrotron emission (which typically has values of 0.7 to 0.8 
for star-forming galaxies (Condon 1992)), and any dependence on 
galaxy properties such as dust temperature; Blain (1999) showed 
that this method would be unable to distinguish hot dusty objects at 
high redshift from cold dusty objects at low redshift. A final prob- 
lem is the possible presence of an AGN in one of the high-redshift 
sources. 

Previous attempts to predict how af^^ should depend on red- 
shift have used samples of low-redshift galaxies that are selected in 
the FIR. Our previous work on the OS SLUGS sample (Vlahakis 
et al. 2005) found a population of cold dusty galaxies (T ~ 20 K) 
which is missing from /i?A5-selected samples, and thus these 'cold' 
galaxies have so far been unrepresented in any analysis of the ef- 
fectiveness of the Carilli & Yun method as a redshift indicator. The 
rest-frame wavelength for any source at moderate redshifts (z < 2) 
detected in the SCUBA surveys will be at a significantly longer 
wavelength than the IRAS bands, so a SCUBA galaxy is likely to 
contain more cold dust than a low-redshift galaxy detected by an 
IRAS survey. 

In this section we use the OS SLUGS galaxies to make a 
new analysis of the af^i" -redshift relation, in a similar way to 
that carried out by Dunne et al. (2000a) for the IRS sample. 
We use the fitted SEDs for the 17 OS galaxies detected at both 
850 [im and 450 vim (fitted using a two-component model to the 
60 \im, 100 |i.m, 450 \im and 850 [im data points as described in 
Vlahakis et al. (2005)). The 1.4 GHz fluxes used are listed in Ta- 
ble [T] The change in af^f with redshift for each of the OS galaxy 
SEDs is plotted in Figure [3] shown as the dashed lines. The thick 
solid lines show the median of the predictions for the 25 IRS galax- 
ies with the highest 60- vim luminosities, together with the ±1(7 
variation in the predictions. The plotted points in the diagram show 
all the deep SCUBA sources with both spectroscopic redshifts and 
radio measurements from Chapman et al. (2005). Stars indicate the 
Chapman et al. sources which have the most robustly determined 
redshifts — following Aretxaga et al. (2007), those whose spec- 
troscopic redshifts are derived from the identification of at least 
two spectral features and which have unambiguous optical/IR/radio 
counterparts. 

It is immediately evident that many of the OS galaxies predict 
a very different afl'l' -redshift relationship to the IRS galaxies, with 
about two thirds of the OS galaxies lying outside the la variation 
of the IRS predictions and having higher af^4 at a given redshift. 
Therefore, in the absence of an independent measurement of red- 
shift, a given source could be either IRS-like and at higher redshift 
or OS-like and at lower redshift, which is essentially the degen- 
eracy noted by Blain (1999). When the OS galaxies are included, 
it can be seen that the uncertainty in the Carilli & Yun method is 
much greater than would be evident from the IRS curves alone. In 
summary, it seems that the temperature of the bulk of the dust sig- 
nificantly affects the position of a galaxy on the a^A-z diagram, 
and that as a consequence of the differing proportions of warm and 
cold dust in local galaxies, which is quite clearly shown in the com- 
parison of the two SLUGS samples, the Carilli & Yun technique 
is unlikely to be reliable as a redshift estimator for deep submm 
sources without measuring the temperature of the dust in these ob- 
jects. 

We now turn to the measurements for the high-redshift galax- 
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Figure 3. Tlie «f -z relationship for the 17 OS SLUGS galaxies detected 
at both 850 and 450 \Mn, shown as the dashed lines. The thick solid lines 
show the median of the predictions for the 25 most luminous IRS sample 
galaxies, together with the itlcr spread for these predictions. The plotted 
points in the diagram show all the deep SCUBA sources with both spectro- 
scopic redshifts and radio measurements from Chapman et al. (2005). Stars 
indicate those which have the most robustly determined redshifts. follow- 
ing Ai'etxaga et al. (2007). The fact that the OS galaxies tend to have higher 
values of a^^^ than the IRS galaxies can clearly be seen, with about two 
thirds of the OS curves having higher af^^ ^ given redshift than the lo" 
variation of the IRS predictions. 

ies. The main thing to note from Figure[3]is that the high-redshift 
sources are significantly brighter at radio wavelengths, or alterna- 
tively fainter at FIR wavelengths, than the predictions for either the 
OS or IRS sample. Lewis, Chapman & Helou (2005) find a sim- 
ilar result comparing their model predictions to measurements of 
deep submm sources. Evidence that the ratio of radio to FIR emis- 
sion may be higher at high redshift than at low redshift has also 
been found by other authors (e.g. Kovacs et al. 2006, who used 
additional submm observations of a subset of the Chapman et al. 
(2005) sample). 

There are a four possible explanations for the difference be- 
tween the Chapman et al. sources and the SLUGS af^i' -z predic- 
tions. Only one is that there is a genuine difference between the 
FIR-radio relationships at high and low redshift. 

(i) The high-redshift SCUBA galaxies have an average dust 
luminosity that is about four times greater than the most lumi- 
nous SLUGS galaxy (Coppin et al., in preparation). Therefore, if 
there is an inverse correlation between af^^ ^nd luminosity at zero 
redshift, our models will under-predict the radio emission of the 
high-redshift galaxies. We investigate this possibility in Figure |4] 
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cause a) the galaxies are not generally strong X-ray sources (Was- 
kett et al. 2003), b) the radio morphologies are not those typical 
of radio-loud AGN (Chapman et al. 2004), and c) the optical spec- 
tra are often those of starbursts rather than AGN (Chapman et al. 
2005). In addition, we find that excluding the 18 submm galaxies 
which exhibit AGN characteristics does not affect the result we see 
in Figure[3] 

(iii) Another possibility, which would contribute to the scatter of 
the Chapman et al. sources in Figure[3l is that some of the redshifts 
are unreliable. In order to test this we selected a subsample of the 
Chapman et al. sources with the most robustly determined redshifts 
— following Aretxaga et al. (2007), those whose spectroscopic red- 
shifts are derived from the identification of at least two spectral fea- 
tures and which have unambiguous optical/IR/radio counterparts. 
We find that using this reduced sample does not affect the result 
we see in Figure [3] and as such it seems unlikely that unreliable 
redshifts are the explanation. 

(iv) The fourth possibility, which we feel is most likely, is that 
the relationship between FIR and radio emission at high redshift 
is different from that at low redshift. Given the very different con- 
ditions in the Universe ten billion years ago compared to today, it 
would actually be quite surprising if the relationship were the same. 
Moreover, there is already evidence that the radio sources associ- 
ated with high-redshift dust sources are quite different from those 
today. Chapman et al. (2004) find that the majority of the high- 
redshift sources have larger radio sizes (~ 10 kpc) than are found 
for local Ultra-Luminous Infrared Galaxies (ULIRGs, < Ikpc). 
This would suggest that the physical conditions may be different 
for low and high redshift sources, in which case one might expect 
to see some evolution of the FIR-radio relation. 



(b) 



Figure 4. Submm-radio spectral index, Qi 4, versus (a) 850- 
luminosity and (b) 60- [am luminosity, for the OS and IRS samples. 
Filled circles indicate OS galaxies detected at 850 |jm while open circles 
((a) only) represent the non-detections at 850 [am; galaxies with both 
submm and radio upper limits are not plotted. The IRS galaxies are plotted 
as cross symbols. 



There is only a weak inverse correlation between af^f and 850- 
^.m luminosity (Figure |4(a)[(. There is an inverse correlation be- 



tween Qi 4 and 60- \im luminosity (Figure |4(b)| (, but it is unclear 
whether it continues at luminosities > lO^'' WHz^^sr^^. The dis- 
crepancy in Figure[3]between our IRS predictions and the Chapman 
et al. sources is about 0.25 in alpha, and the discrepancy is greater 
for the OS predictions. In order for even the smaller of the discrep- 
ancies to be explained by the correlation in Figure [4(b)| ot^2 would 
have to decrease by approximately 0.25 over a luminosity range of 
a factor of ~4. This seems unlikely looking at the slow change in 
a, but we caimot exclude the possibility. 

(ii) One possibility could be that the Chapman et al. sources 
have significant AGN activity contributing to their 1.4 GHz emis- 
sion (an extreme example in the local Universe would be NGC 6166 
and the other six or so outliers in the OS sample FIR~radio correla- 
tion shown in Figure [T]( — if dust is heated to higher temperatures 
in stronger radiation fields, e.g. near to an AGN, we could expect 
to see both higher 1.4 GHz emission and higher dust temperature. 
However, it seems unlikely that the radio and FIR emission from 
these galaxies is dominated by radiation from an active nucleus be- 



4 DISCUSSION 

Although the explanation of the relation between radio and FIR 
emission at low redshift is still uncertain (see SectionfTJ, there are 
still some general statements one can make about the possibility of 
this relationship evolving with redshift. 

We will assume that the standard model is correct: the FIR 
emission is from dust heated by young OB stars and the radio 
emission is from relativistic electrons which were produced by su- 
pernovae. Within the standard model (which invokes an explicit 
dependence on star formation), there are two explanations of the 
slope (of roughly unity on a log-log plot) of the low-redshift re- 
lationship. The first is the calorimeter model (Volk 1989). In this 
model, the electrons radiate all their energy before they escape from 
the galaxy; variations in the magnetic field change the rate at which 
an individual electron loses its energy but do not change the FIR- 
radio ratio. The second type of explanation is that the slope, which 
is not precisely unity, is a cosmic conspiracy. In models of this type, 
electrons do escape from galaxies before they radiate all their en- 
ergy, which would lead to a lower FIR-radio ratio in galaxies with 
stronger magnetic fields. However, this is avoided in the models 
by fine-tuning them so that, for example, electrons diffuse more 
rapidly out of galaxies in which the magnetic fields are stronger 
(e.g. Chi & Wolfendale 1990). Other fine-tuning models include, 
for example, those of Helou & Bicay (1993) and Niklas & Beck 
(1997). The calorimeter model seems the more natural explanation, 
but its disadvantage has been that the radio spectra of galaxies do 
not show the clear down-turn at high frequencies expected for an 
aging population, although Thompson et al. (2006) have argued that 
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this is the result of another cosmic conspiracy - of ionization and 
bremsstrahlung losses flattening the spectra. 

If the calorimeter model is correct, we would probably expect 
a change in the FIR-radio ratio with redshift in the opposite direc- 
tion to the one we actually observe. In this model, increases in the 
magnetic field in high-redshift galaxies, which one might expect as 
the result of increased turbulence in the interstellar medium, will 
have no affect on the FIR-radio ratio. The one effect, however, that 
one might expect is that the electrons will lose a greater fraction 
of their energy to inverse Compton losses, because of the increase 
in the energy density of the microwave background and because of 
the strong interstellar radiation field in these galaxies. In this case, 
however, we would expect the FIR-radio ratio to increase with red- 
shift because an electron is losing more of its energy by inverse 
Compton losses than by synchrotron losses. 

The cosmic conspiracy models contain more free parameters 
and therefore lead to less definite predictions. However, since the 
high-redshift galaxies detected in the deep submm surveys do ap- 
pear to have larger radio sizes than the ULIRGs found in the local 
Universe, one might expect the relativistic electrons to escape more 
slowly from the galaxies. In this case, one would expect the FIR- 
radio ratio to decrease with redshift. This is in qualitative agree- 
ment with what we observe. Therefore, our results provide some 
evidence against the calorimeter model. 

One important question is whether the SLUGS SEDs are suit- 
able templates for the high-redshift objects in the Chapman et al. 
(2005) sample. Unfortunately, for these distant submm sources 
there are very limited observations at wavelengths which probe the 
peak in the dust SED and the spectral slope in the radio regime. 
The current best estimates of the FIR SEDs come from SHARC-2 
observations at 350nm (probing rest frame 120vim at z = 2) pre- 
sented by Kovacs et al. (2006) and Coppin et al. (in preparation). 
This allows a basic estimate of the dust temperature independent of 
assumptions about the FIR-radio correlation, and these authors find 
an average dust temperature of ~ 35 K for submm selected galax- 
ies — very similar to the average temperature for the SLUGS IRS 
sample (Dunne et al. 2000b). In this respect, the FIR/submm part of 
the IRS SLUGS SEDs are probably a sensible choice of template. 
In the radio, however, it is less obvious that things could be the 
same at both low and high redshift. For example, there are very few 
measurements of the radio spectral index for high-redshift submm 
sources. A flattening of the average spectral index compared to that 
for the low-redshift sources (a = 0.7) would tend to move the 
high-redshift sources in the direction observed in Figure [3] Hunt, 
Bianchi & Maiolino (2005) notice that SEDs appropriate for com- 
pact dwarf galaxies (which have a higher component of free-free 
emission and so flatter spectral indices) produce better photometric 
redshift estimators than SEDs based on more typical IR luminous 
galaxies. The reason why the high-redshift sources are generally 
brighter at radio wavelengths compared to our SLUGS af ^4 -z pre- 
dictions may be uncovered by deep multi-wavelength radio obser- 
vations of submm sources which will investigate their radio proper- 
ties, as well as better coverage of the IR SED (e.g. from Herschel). 



5 CONCLUSIONS 

Using the results of the OS SLUGS sample we have investigated the 
relationships between FIR-submm and radio properties for 'nor- 
mal' local galaxies, compared to bright IRAS galaxies (using the 
IRS SLUGS sample from Dunne et al. 2000b). The results for the 
OS galaxies have then been used to assess the reliability of the Car- 



illi & Yun (1999, 2000) radio-submm redshift estimator technique, 
for the first time using a sample containing significant fractions of 
cold dust. 

We find the following results. 

• There is a very tight correlation between 1.4 GHz luminos- 
ity and 60-|J.m luminosity for the OS sample galaxies. There 
is much more scatter in the correlation between radio and 
850 \im luminosity, thus providing evidence that the standard ex- 
planation for the FIR-radio luminosity correlation is correct and 
that it is the population of young (OB) stars which both heats the 
dust and provides the source of radio emission. 

• Using the SEDs of the 17 OS galaxies detected at both 850- 
and 450- [im we have shown that for a sample of 'normal' galaxies 
there is much more scatter in the af^f -redshift relation than seen 
for the 25 most luminous IRAS galaxies of the IRS SLUGS sam- 
ple, most likely since the OS objects contain large fractions of cold 
dust. From this we conclude that in order for the Carilli & Yun 
radio-submm redshift indicator method to be reliable as a redshift 
estimator for deep submm sources one would need measurements 
of the temperature of the dust in those objects. 

• We have compared our OS af ^4 -redshift relations to the deep 
submm observations of Chapman et al. (2005) and demonstrate that 
the large majority of these deep SCUBA sources must have very 
different properties to our sample of low-redshift SLUGS galaxies. 
We find that the majority of the deep SCUBA galaxies are brighter 
at radio wavelengths than any of our predictions. We cannot ex- 
clude the possibility that this simply reflects that our low-redshift 
galaxies are not suitable templates for the high-redshift SCUBA 
galaxies, but we think it more likely that this shows the FIR-radio 
relation is evolving with redshift. 
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